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ABSTRACT 


The  pole  face  losses  associated  with  the  alternating  flux 
density  that  is  produced  by  the  slot  and  teeth  in  the  armature 
for  the  467  kva  experimental  generator  were  experimentally 
evaluated  by  using  an  eddy-current  dynamometer.  The  design 
factors  which  influence  pole  face  losses  in  the  generator  were 
duplicated  in  the  test  fixture  as  close  as.  possible.  Centrif¬ 
ugal  stress  in  the  rotor  limited  the  tooth  ripple  frequency 
to  one-half  the  tooth  ripple  frequency  of  the  generator  and 
it  was  necessary  to  extrapolate  the  data  from  24,000  cps  to 
48,000  cps  to  determine  the  pole  face  losses  at  these  conditions. 
Ho  attempt  was  made  to  control  the  temperature  of  the  pole  faces. 
Therefore,  the  test  results  had  to  be  corrected  for  a  change  of 
resistivity  since  the  pole  faces  will  operate  at  a  much  higher 
temperature  in  the  generator. 

A  complete  analysis  of  the  pole  face  loss  problem  based 
upon  theoretical  considerations  is  presented  in  the  report. 

The  problem  of  dealing  quantitatively  with  the  magnetic  flux 
density  distribution  at  the  surface  of  the  poles  is  treated 
analytically.  Harmonic  frequencies  and  the  losses  due  to 
each  harmonic  frequency  were  considered.  The  effect  of  grooves 
in  the  pole  faces  were  also  considered  in  the  theoretical 
analysis.  The  theoretical  equations  show  that  the  eddy-current 
losses  are  a  function  of  permeability  of  the  pole  face  material. 
In  the  actual  problem  this  quantity  is  a  variable  and  a  method 
for  determining  the  applicable  value  of  this  variable  is  pre¬ 
sented  . 

The  general  conclusions  derived  from  this  work  are:  (1) 
materials  which  offer  the  best  combination  of  d~c  magnetic  and 
strength  properties  for  a  high  speed,  high  temperature,  solid 
rotor  have  inferior  properties  from  the  standpoint  of  pole 
face  losses  and  (2)  a  16-fold  decrease  of  eddy-current  losses 
can  be  obtained  by  using  a  laminated  material  for  the  pole 
tips  which  has  been  chosen  for  its  high  a-c  magnetic  permeability 
when  compared  with  the  losses  associated  with  material  that  has 
been  selected  for  its  creep  strength  or  other  criteria  except 
its  a-c  magnetic  permeability. 

{This  abstract  is  subject  to  special  export  controls  and 
each  transmittal  to  foreign  governments  or  foreign  nationals 
may  be  made  only  with  prior  approval  of  the  Air  Force  Aero 
Propulsion  Laboratory,  APIP-1.) 
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SECTION  I 


INTRODUCTION 


Accurate  evaluation  of  the  power  losses  in  the  rotor  pole 
faces  is  required  to  determine  the  thermal  state  of  a  rotor. 

The  -mechanical  integrity  of  the  rotating  mass  depends  upon 
these  losses  since  the  allowable  stresses  depends  upon  the 
thermal  state  of  the  rotor.  Losses  which  produce  heat  energy 
in  the  pole  faces  are  caused  by  two  sources  of  flux  oscilla¬ 
tions.  The  first  is  due  to  the  fact  that  as  the  pole  face 
passes  under  a  slotted  armature  all  points  on  the  pole  face 
experience  a  varying  flux  density.  At  one  instant  a  point 
on  the  pole  face  is  under  a  tooth  where  the  radial  component 
of  the  flux  density  is  a  maximum  and  at  the  next  instant  it 
is  under  a  slot  where  the  radial  component  of  the  flux  density 
is  a  minumum.  The  varying  flux  density  in  the  radial  direction 
within  the  poles  is  associated  with  a  varying  flux  in  the 
tangential  direction  since  all  flux  lines  must  form  a  closed 
loop.  The  varying  flux  density  produces  an  emf  which  gives 
rise  to  eddy-currentc  that  produce  losses  in  the  pole.  The 
varying  flux  also  produces  minor  hysteresis  loops  which  cause 
additional  losses.  The  losses  produced  by  the  slots  and  teeth 
in  the  armature  are  generally  called  "tooth  ripple  losses". 

These  losses  exist  whenever  the  machine  is  magnetized. 

A  second  source  of  losses  in  the  poles  is  produced  by  the 
varying  flux  that  is  generated  by  harmonic  mmf ’s  associated 
with  the  armature  windings.  These  losses  exist  only  when  current 
is  flowing  in  the  windings. 

Many  investigators  have  published  equations  for  evaluating 
the  power  losses  due  to  eddy-currents  and,  for  the  most  part, 
they  have  been  collated  in  a  paper  by  Aston  tind  Rao  (reference 
1).  Equations  for  the  evaluation  of  the  losses  associated  with 
the  varying  flux  produced  by  the  armature  windings  are  presented 
in  Reference  2.  Most  investigators  have  used  semi-empirical 
or  purely  empirical  means  based  upon  the  results  of  several 
experiments  to  establish  their  formulas.  Other  investigators 
have  based  their  equations  upon  theoretical  considerations 
by  admission  of  certain  approximations.  In  general,  the  theore¬ 
tical  treatments  assume  that  only  eddy-current  loss  has  impor¬ 
tance  and  is  always  predominant  over  hysteretic  loss.  Admission 
of  approximations  are  necessary  to  deal  with  the  permeability 
of  the  material,  which  is  a  variable  quantity  and  is  always 
treated  as  a  constant.-  Admission  of  certain  approximations 
are  also  necessary  to  obtain  numerical  results  for  the  distri¬ 
bution  of  flux  density  on  the  surface  of  the  pole  faces. 
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The  value  of  the  parameters  which  influence  pole  face 
losses  for  the  467  kva  experimental  generator  is  far  above 
the  range  previously  investigated.  Hence,  empirical  equations 
which  contain  numerical  constants  that  are  based  upon  experi¬ 
mental  results  from  average  designs  are  apt  to  give  erroneous 
results  for  the  novel  467  kva  generator.  Hence,  equations 
which  are  based  non  theoretical  considerations  are  preferred, 
since  this  method  can  be  employed  for  a  wider  range  of  values 
of  the  variables.  One  of  the  important  variables  which  has  a 
significant  influence  upon  eddy-current  losses  is  the  resis¬ 
tivity  of  material  in  the  poles.  The  pole  fa  ,es  in  the  experi¬ 
mental  generator  will  operate  at  approximately  10C0°F  and  the 
material  in  the  poles  will  have  a  higher  resistivity  than  the 
average  design.  Therefore,  experimental  results  obtained  with 
average  designs  can  not  be  used  to  predict  the  eddy-current 
losses  in  a  high  temperature  generator  unless  the  formulas 
correctly  describe  the  influence  of  resistivity.  The  experi¬ 
mental  generator  tooth  ripple  frequency  will  be  48,000  cps, 
which  is  approximately  eight  times  greater  than  the  frequency 
previously  investigated.  Thus,  formulas  based  upon  experi¬ 
mental  results  from  low  frequency  tests  must  correctly  describe 
the  influence  of  frequency  before  they  can  be  used  to  predict  * 

the  losses  for  the  experimental  generator. 

Several  of  the  available  equatipns  were  used  to  evaluate  » 

the  tooth  ripple  losses  for  the  generator  and  the  results  were 
presented  in  Reference  3.  A  wide  disagreement  was  found  to 
exist . 

In  general,  pole  face  losses  are  of  minor  importance  in 
slow-speed  machines  and  an  order  of  magnitude  difference  be¬ 
tween  .the  calculated  and  actual  values  may  be  cf  little  signifi¬ 
cance  with  regard  to  machine  performance.  For  the  high  speed 
experimental  generator,  where  the  rotor  material  is  being  used 
near  its  maximum  stress  and  temperature  capabilities,  an  error 
of  less  than  50  percent  could  result  in  a  complete  failure  of 
the  machine.  For  this  reason  and  since  the  accuracy  of  avail¬ 
able  equations  is  questionable,  it  was  deemed  necessary  to 
perform  experiments  to  determine  the  tooth  ripple  losses  with 
material  and  design  factors  as  close  as  possible  to  those 
contemplated  for  the  generator.  The  design  of  a  pole  face 
loss  test  fixture  was  initiated  in  1963.  The  fixture  was  built 
and  tests  were  accomplished  during  1964  and  reported  in  Refer¬ 
ence  3.  In  general,  the  results  from  these  tests  were  incon¬ 
clusive  with  regard  to  pole  face  losses.  The  test  fixture  and 
instrumentation  was  modified  to  reduce  the  experimental  error 
and  the  program  was  continued.  In  addition,  a  complete  theore¬ 
tical  analysis  of  the  problem  was  performed.  The  results  from 
this  analysis  and  the  experimental  effort  are  presented  in 
in  this  report. 
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SECTION  II 


TEST  FIXTURE 


The  distinguishing  features  of  the  test  fixture  employed 
for  this  investigation  are  illustrated  by  Figure  1.  In  effect, 
the  -fixture  was  an  eddy-current  dynamometer  or  a  magnetic 
brake  in  which  a  toothed  rotor  was  rotated  between  two  poles 
of  an  electromagnet.  Torque  reaction  on  the  trunnion  mounted 
stator  was  measured  on  a  scale  to  obtain  the  electrical  and 
windage  losses  in  the  machinery.  The  air  gap  was  controlled 
by  adjusting  the  radial  position  of  the  poles. 

The  tooth  pitch  and  slot  configuration  in  the  rotor  dupli¬ 
cated  the  tooth  pitch  and  slot  configuration  of  the  armature 
for  the  467  kva  experimental  generator  shown  in  Figure  2. 
Centrifugal  stress  limited  the  speed  of  the  120- tooth  rotor 
to  12,000  rpm.  Thus,  the  maximum  tooth  ripple  frequency  that 
could  be  attained  was  24,000  cps  or  one-half  the  tooth  ripple 
frequency  of  the  generator. 

The  pole  specimens  used  in  the  fixture  are  shown  in  Figures 
3a  and  3b.  The  specimens  were  made  of  E-ll  steel  and  were  heat 
treated  to  hardness  specified  for  the  full  scale  rotor.  The 
grooves  in  the  pole  faces  duplicated  the  grooves  used  in  the 
rotor  pole  faces  except  for  width;  they  were  approximately 
0.010  inches  wide  at  the  surface  of  the  pole  or  approximately 
0.005  inches  wider  than  desired. 

It  would  have  been  highly  desirable  to  construct  a  test 
fixture  in  which  no  losses  except  the  tooth  ripple  losses 
could  have  existed.  In  general,  this  is  not  capable  of  being 
realized.  Designers  of  a  pole  face  loss  test  fixture  must  cope 
with  the  problem  of  separating  the  pole  face  losses  from  the 
other  losses  in  the  machine.  Ingenious  methods  are  always 
required.  The  losses  in  the  test  fixture  shown  in  Figure  1 
which  caused  a  torque  reaction  on  the  scale  were:  (1J  the 
loss  associated  with  the  aerodynamic  drag  on  the  stator, 

(2)  rotor  iron  losses,  (3)  magnetic  circuit  pulsation  loss, 
and  (4)  pole  face  losses.  The  reaction  in  the  bearings  which 
support  the  high  speed  rotor  and  the  aerodynamic  reaction  on 
the  rotor  were  not  imposed  upon  the  trunnion  mounted  stator. 

The  initial  tests  performed  in  1964  (reference  3)  showed  that 
the  aerodynamic  reaction  on  the  stator  was  approximately  ten 
times  greater  than  the  reaction  due  to  electrical  losses  in 
the  machine  at  operating  conditions  pertinent  to  the  generator. 
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Figure  1.  Pole  Face  Loss  Test  Fixture 
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Stationary  shrouds  were  installed  on  the  windward  side  of  the 
poles  to  decrease  the  aerodynamic  reaction  on  the  stator  after 
the  initial  tests  were  performed.  With  this  modification  the 
windage  reaction  was  decreased  to  approximately  two  times  the 
reaction  due  to  electrical  losses  at  conditions  relevant  to  the 
actual  generator.  The  windage  reaction  was  measured  by  noting 
the  scale  reading  with  the  fixture  demagnetized  for  each  test 
speed  and  air  gap.  This  value  was  then  subtracted  from  the 
readings  taken  with  the  fixture  magnetized.  The  net  reaction 
represented  the  electrical  losses  in  the  magnetic  circuit. 
Magnetic  circuit  pulsation  losses  could  not  be  separated  from 
the  other  losses.  Pulsation  losses  are  caused  by  a  variation 
of  air  gap  reluctance  and  are  theoretically  present  when  the 
pole  arc  covers  a  non-integral  number  of  slot  pitches  and  is 
most  severe  when  the  fractional  part  of  the  number  is  one-half 
for  an  integral  slot  winding.  Since  the  test  fixture  had  an 
armature  with  an  integral  slot  winding  feature  or  sixty  slot 
pitches  per  pole  pitch,  the  number  of  slot  pitches  per  pole  arc 
was  made  an  integral  number,  or  fifteen  in  order  to  reduce  the 
pulsation  losses. 

To  verify  that,  the  test  fixture  design  was  such  that  pulsa¬ 
tion  losses  were  negligible,  a  forty  turn  coil  was  wound  on  the 
stator  frame  so  that  any  flux  pulsation  in  the  magnetic  circuit 
could  be  measured.  This  was  accomplished  by  connecting  a  cali¬ 
brated  oscilloscope  to  the  coil  and  measuring  the  induced  voltages 
and  frequency.  Figure  4  shows  the  induced  coil  voltage  as  taken 
from  polaroid  pictures;  curves  are  presented  for  several  air  gaps. 

The  maximum  magnitude  of  the  flux  pulsation  for  the  smallest 
air  gap  was  approximately  2  lines  per  square  inch  and  it  is 
reasonable  to  conclude  that  the  loss  produced  by  this  magnitude 
of  pulsation  was  very  minute. 

The  electrical  losses  in  the  fixture  were  separated  into 
rotor  iron  loss  and  pole  face  loss  by  determining  the  rotor 
iron  loss  for  a  given  rotor  speed  and  frame  flux  at  large  air 
gaps.  Since  pole  face  loss  is  a  function  of  air  gap,  this 
loss  decreased  when  the  air  gap  increased. 

Since  rotor  iron  loss  is  independent  of  air  gap,  it  was 
obtained  by  plotting  the  measured  electrical  losses  for  constant 
values  of  flux  and  speed  as  a  function  of  air  gap  and  extrapola¬ 
ting  the  curves  to  infinite  air  gap  where  pole  face  losses 
are  zero.  The  rotor  iron  loss  was  then  subtracted  from  the 
total  electrical  loss  to  obtain  the  pole  face  loss  for  a  given 
frame  flux,  air  gap,  and  rotor  speed. 
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SECTION  III 


TEST  RESULTS 


The  initial  test  runs  had  indicated  that  a  considerable 
amount  of  flux  was  fringing  around  the  pole  face  flux  coils 
and  that  the  total  flux  through  the  rotor  was  higher  than  that 
indicated  by  the  pole  face  coil  readings.  To  verify  this, 
another  set  of  coils  were  wound  on  the  stator  frame  so  that 
all  of  the  flux  entering  the  rotor  could  be  measured.  Figure 
5  shows  a  plot  of  measured  flux  as  a  function  of  the  field  coil 
excitation  current  for  both  the  pole  face  and  frame  flux  coils 
for  the  case  of  a  radial  gap  of  0.092  inches.  These  curves 
verify  the  conclusion  that  the  total  rotor  flux  is  significantly 
higher  than  the  flux  in  the  stator  pole  faces. 

Curves  of  rotor  flux  and  pole  face  flux  plotted  as  a 
function  of  field  coil  current  for  the  four  test  air  gap  are 
shown  in  Figures  6  and  7  respectively.  These  curves  are  inde¬ 
pendent  of  rotor  speed. 

One  of  the  initial  test  runs  included  an  air  gap  of  0.054 
inch  but  it  was  later  found  when  plotting  the  data  that  an 
insufficient  number  of  data  points  were  taken  for  rotor  flux 
values  between  0  and  1000  kilolines  (only  two  points  were 
included).  Because  of  the  steep  slope  of  the  air  gap  line 
for  this  small  gap  length,  smaller  increments  of  field  current 
should  have  been  used  to  obtain  more  points.  For  this  reason, 
test  data  for  the  0.054  inch  air  gap  was  omitted  from  the 
analysis. 

Figures  8,  9,  and  10  show  plots  of  measured  test  fixture 
electrical  loss  (i.e.,  total  measured  loss  minus  windage  loss) 
plotted  as  a  function  of  field  coil  current  for  the  rotor  speeds 
of  4,000  rpm,  8,000  rpm,  and  12,000  rpm  respectively.  The  high 
speed  limit  for  the  test  fixture  is  12,000  rpm. 

From  the  curves  of  rotor  flux  versus  field  current  and 
electrical  loss  versus  field  current  for  the  various  test  air 
gaps,  plots  of  electrical  loss  as  a  function  of  air  gap  length 
for  constant  values  of  rotor  flux  were  obtained.  These  curves 
are  shown  in  Figures  11,  12,  and  13  for  rotor  speeds  of  4,000 
rpm,  8,000  rpm,  and  12,000  rpm  respectively.  The  rise  in 
electrical  loss  with  smaller  air  gaps  due  to  increased  pole 
face  loss,  is  very  evident.  The  right  hand  side  of  the  curves, 
where  the  rate  of  change  of  loss  with  air  gap  is  practically 
zero,  represents  the  rotor  tooth  and  core  loss.  The  difference 
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between  the  loss  values  given  by  the  curves  at  a  particular  air 
gap  and  the  rotor  tooth  and  core  loss  represents  the  fixture 
pole  face  loss.  Pole  face  loss  is  a  function  of  flux  in  the 
pole  races  and  not  the  total  flux  in  rotor.  The  flux  in  the 
pole  faces  for  a  given  value  of  rotor  flux  was  determined  from 
the  curves  of  Figures  6  and  7.  Pole  face  loss,  expressed  in 
units  of  watts  per  square  inch,  and  plotted  as  a  function  of 
pole  face  flux  density  ana  tooth  frequency,  are  shown  in  Figures 
14,  15.  and  15.  The  units  of  flux  density,  and  tooth  frequency 
were  used  since  these  units  are  more  descriptive  of  the  factors 
involved  in  the  creation  of  pole  face  loss.  21  tooth  frequency 
of  24,000  cps  corresponds,  in  this  case,  to  a  rotor  speed  of 
12,000  rpm. 

Effects  of  tooth  saturation  are  clearly  evident  on  the  test 
curves  as  sncwn  by  the  decreasing  slope  of  the  loss  curve  for 
pole  face  flux  densities  above  approximately  60  kilolines  per 
square  inch.  Saturation  of  the  teeth  decreases  the  amplitude 
of  the  flux  ripple  in  the  pole  surface  thereby  reducing  losses. 

The  radial  gap  length  for  the  experimental  generator  is 
0.100  inch.  Since  test  data  was  not  obtained  for  this  exact  gap 
length,  the  loss  corresponding  to  a  0.100  inch  air  gap  was 
obtained  by  plotting  the  test  loss  values  versus  the  air  gap 
lengths  tested  as  shewn  in  Figure  17.  Figure  IS  shows  a  plot 
of  pole  face  loss  as  a  function  of  tooth  frequency .  Loss  values 
derived  from  the  test  data  axe  shown  by  the  circled  points.  The 
tooth  frequency  for  the  experimental  generator  is  48 . 000  cps. 

It  was  therefore  necessary  to  extrapolate  these  data  points  to 
48,000  cps  as  indicated  in  Figure  IS.  The  accuracy  of  the 
extrapolated  loss  is  sensitive  to  errors  in  the  data  points 
that  were  obtained  from  the  experimental  values.  Since  the 
torque  on  the  trunnion  mounted  stator  due  to  pole  face  loss 
was  small  in  comparison  with  the  combined  torque,  the  accuracy 
of  torque  measurement  had  a  predominant  influence  upon  accuracy 
of  the  experimental  data.  Measurement  error  associated  with 
the  torque  on  the  stator  was  approximately  +  one  percent.  For  a 
pole  face  flux  of  56.5  kilolines  per  square  inch  the  windage  loss 
torque  was  less  than  or  equal  to  200  percent  of  the  electrical 
loss.  Thus  the  maximum  error  associated  with  the  electrical  loss 
which  was  obtained  by  substracting  the  windage  loss  was 

erior  =  {1+0. Qi)  Cl+W/Ej  -  {i~G.0X}R/E  -I  =  0.02W/S+0.01 

where  W  =  windage  loss,  E  =  electrical  loss.  For  W/E  =  2.0  the 
error  becomes  +  five  percent.  The  error  associated  with  the  pole 
face  loss,  which  was  obtained  by  subtracting  the  electrical  loss 
for  a  large  air  gap,  becomes 

error  =  (1+G.05)E/P  -  (I~0,05jS/P  -  1  --  0.1S/P  +  0.05 
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where  P  =  pole  face  loss  and  E  =  electrical  loss  for  a  large  air 
gap.  For  an  air  gap  of  approximately  0.1  inch  the  value  of  E/P 
was  about  two.  Therefore,  the  error  associated  with  the  pole 
face  losses  represented  by  the  three  points  in  Figure  18  could 
have  been  +  25  percent.  The  final  experimental  points  shown 
in  Figure  T8  are  probably  more  accurate  since  conformity  of  the 
experimental  values  was  obtained  by  plotting  various  dependent 
quantities  against  the  independent  variables  of  the  experiment. 

In  general,  this  allowed  values  for  the  dependent  variable  to  be 
obtained  from  curves  which  represented  the  locus  of  the  experi¬ 
mental  values  and  eliminated  part  of  the  error  caused  by  measure¬ 
ment  accuracy.  Unfortunately,  the  accuracy  of  the  windage  torque 
data  could  not  be  improved  since  windage  torque  measurements 
could  only  be  obtained  before  and  after  the  test  fixture  was 
magnetized  for  a  particular  calibration.  These  readings  varied 
by  +  one  percent.  When  the  combined  torque  values  were  taken 
from  curves  without  the  normal  variation  of  points,  the  accuracy 
of  the  combined  torque  data  was  improved  from  +  one  percent  to 
+0.2  percent.  Therefore,  the  error  associated  with  the  electri¬ 
cal  loss  was  reduced  from  +  five  percent  to  +  2.6  percent  for 
cases  with  a  windage  loss  Tess  than  or  equal  to  200  percent  of 
the  electrical  loss.  The  error  associated  with  the  pole  face 
loss  was  therefore  reduced  from  +  25  percent  to  +  13  percent 
for  cases  with  a  small  air  gap.  Nevertheless,  tEe  possibility 
of  a  larger  error  exists  fcr  the  point  given  for  48,000  cps  in 
Figure  18  since  the  loss  curve  had  to  be  extrapolated  from 
24,000  to  48,000  cps  using  three  data  points  which  may  be  in 
error  by  13  percent.  Extrapolations  on  a  logarithmic  plot  were 
made  using  the  three  data  points  given  in  Figure  18  with  +  13 
percent  variation  of  the  values.  A  straight  line  extension 
through  the  maximum  point  at  24,000  cps  and  the  minimum  point 
at  8,000  cps  and  vice  versa  resulted  in  a  variation  of  32  per¬ 
cent  at  48,000  cps.  It  is  logical  to  assume  that  this  variation 
represents  the  possible  error  for  the  extrapolated  value  at 
48,000  cps.  The  extrapolation  gives  a  pole  face  loss  of  31.7 
+  10.1  watts  per  square  inch  for  a  pole  face  temperature  of  75°F. 
A  temperature  calculation  performed  for  the  experimental  genera¬ 
tor  with  an  average  coolant  temperature  of  600°F  gave  a  pole  face 
temperature  of  about  980°F.  The  loss  decreases  with  increasing 
temperature  due  to  a  change  in  the  resistivity  of  the  pole  face 
material.  A  correction  for  a  change  in  resistance  as  outlined 
in  Section  V  gives  23.0  +  7.4  watts  per  square  inch  for  the  pole 
face  loss  at  980°F.  The  pole  face  area  for  the  experimental 
generator  is  77.5  square  inches  and  when  multiplied  by  the 
above  loss  density  results  in  a  total  loss  of  1783  +  590  watts. 
This  is  the  loss  component  induced  in  the  rotor  pole  faces  by 
the  flux  "bundling"  effect  of  the  armature  teeth. 


26 


Figure  18.  Pole  Face  Loss  Versus  Tooth  Ripple  Frequency  for 
Radial  Gap  of  0.100  Inch  and  a  Pole  Face  Flux 
Density  of  56.5  Kilolines  per  Square  Inch 
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SECTION  IV 


THEORETICAL  TREATMENT  OF  POLE  FACE  LOSSES 


Eddy-Current  Loss 


-  In  practice  the  characteristic  dimensions  of  the  slots  as 
shown  by  Figure  19a  are  small  when  compared  with  the  dimensions 
of  the  rotor.  Therefore,  a  mathematical  treatment  of  the  pro¬ 
blem  shall  be  made  using  rectangular  coordinates  as  illustrated 
in  Figure  19b. 


The  field  equations  can  be  used  to  analytically  describe 
the  flux,  current,  and  the  eddy-current  loss  in  the  pole  mater¬ 
ial.  Maxwell's  equation  from  Faraday's  law  gives: 

3B 

7XE  =  -ft 

and  the  differential  form  of  Maxwell's  equation  derived  from 
Ampere's  law  is: 

JD 

VXH  »  4*  (J  +  "at ) 

Maxwell’s  magnetic  field  equation  as  derived  from  Gauss's  law 
gives  the  differential  or  point  relation: 


7*3  «  0 


Thus  far,  we  have  stated  three  of  Maxwell's  four  equations  which 
apply  at  a  point  in  a  time- changing  field.  The  fourth  is  Max¬ 
well's  electric  field  equation  as  derived  from  Gauss's  law: 


7*D  =  Charge  density 

Since  we  are  concerned  with  the  eddy-current  loss  in  a  conduct¬ 
ing  medium  we  shall  assume  that  the  medium  is  free  of  charge 
and  has  zero  dielectric  constant.  Thus,  the  displacement 
current : 

,  .  dD  dE 

J  disp.  =  £  at  =  0 

With  this  assumption  the  field  equations  in  electromagnetic 
units  are: 

_  3  B  r.  , 

4  it  J=  7  XH  (2) 

7*  B  =  0  (3) 
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Along  with  Maxwell's  equations  certain  other  fundamental  rela¬ 
tions  are  required  to  establish  the  eddy-current  at  a  point 
within  the  system.  They  are: 


(15) 


Equation  (6)  shows  that  V»J 

V2J 


0  .  Therefore 
3 J 

""p"1”  3t 


writing  Equation  (15)  in  terms  of  the  three  rectangular  components 
of  J  gives  an  equation  which  is  the  vector  sum  of  three  scalar 
equations.  Hence, 

32  Jy  +  32Jy  +  3  2  Jx  =  llIL  3Jy  (16) 

3  x*  3  y~  3  7?  p  ’  3t 

32  Jy  4.  32 Jy  4.  3 2 J y  =  4frp  3 J y  (17) 

r?  3  y  2  3  z2  p*  3t 

32JZ  +  32JZ  4  32 J z  -  9 z  (i8^ 

3  x2  3  y2  3  zr  p  ’  3t 

a  set  of  scalar  equations  can  he  developed  from  Maxwell’s  first 
curl  equation  which  gives  a  relationship  between  the  rectangular 
components  of  B  and  J.  Since, 


VXE  =  vxp’jr-i8  (I9) 

3t 

we  can  express  J  and  B  in  terms  of  the  rectangular  components  to 
give: 


_1 
p 1 

ft  CiBx  +  j By  +  kBz)  = 

J  z 
x(^y 

£*>♦ 

4,3  J  X  _ 
j(9z 

3J  Z-J4 

3X  1 

- 

Ms)* 

(20) 

3X 

3y 

Equation  (20) 

is  the  vector  sum 

of  three 

scalar 

equations. 

Hence, 

1 

3  8  y 

r  z  _  3J  y 

(21) 

"  p* 

3t 

3y  3z 

1 

3  B  y 

=  3J  x  _  3 J  z 

(22) 

F*~ 

3t 

3z  3x 

(23) 

1 

-  P- 

3B  z 
3t 

3J  y  3  J  X 

3x  3y 

In  addition  we  have  from  Equation  (6) 


V  •  J  =  0 


3  J  y  3  J  y  3  J 
3x  +  3y  3z 


(24) 


Equations  (11)  through  (13)  are  differential  equations  for  B 
in  terms  of  the  independent  variables  (x,  y,  z,  and  t)  of  the 
system.  A  solution  for  the  distribution  of  current  density 
must  start  with  a  solution  of  equations  (11)  through  (13) 
since  we  can  define  the  boundary  values  for  the  flux.  Once 
we  have  developed  the  flux  distribution,  we  can  then  develop 
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the  current  distribution  by  using  equations  {16}  through  (18/ 
and.  equations  {21}  through  (24). 


and 


-a2  .  a2  -  iiE  a3  ,  0 

Equations  (27) ,  (28)  ,  and  (29)  are  ordinary  differential  equa¬ 
tions.  In  our  case  the  constants,  a^,  a-2'  a3  ecrual 

=  -a2  =  4tt2  ro2/X2  (30) 

=  a,=  inm2uf  (31) 

T  3 

~  =  a ?=  4ir2m2/X2  +  i8ir2nm}jf/p  *  (32) 


The  solutions  of  these  equations  for  the  case  where  n  =  1  sure: 


X  =  Axexp(i2irmx/X)+Biexp(-i2Trinx/X) 

Y  -  Ciexp^rm/taJy/Xj+D^expC-STfaj''^  y/X) 
T  -  E1exp(im2irft) 


where 


“1 


=  a  + 


P% 


The  solutions  of  Equations  (27) 
where  n  =  -1  are 


(28) ,  and  (29)  for  the  case 


X  =  A_^exp(i2irmx/X)+B_^exp(-i2Tris3c/,\) 

Y  =  C_jexp  (2  wn/w  ^y/x )  -fD_^exp  ( -2irm»^oZJy/X ) 
T  =  E_^exp  (-i27i  left) 


where 

oi_^  =  l-ip2f  X2/p '  »  conjugate  of  =  caj 
For  the  case  where  n  =  1,  we  have 

B*ym  =  [A^expC^Ttax/xl+B^expC-^Tunx/X)]  [C1exp(2rav^y/X)  + 

Biexp(-27Tm>/u2,y/X) J  Exexp(im2^ft)  (33) 

and  for  the  case  whe  *e  n  =  -1,  we  have 

B*vm  =  lA-2exp(i27Tiax/x)+B«1exp(-i27rmx/x)  3  tC-.iexp(2-frm/SiV/x)  + 

D_]^exp (-27rm^r_y/x)  3  E«iexp<-ia2irft)  (34) 
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We  must  now  examine  the  complex  form  of  and 
In  order  for  -+Q  when  y  -*»,  it  is  necessary  that 

exp(*4iiy)  =  expf  (-Um±iVm)y]  and  exp(*/Biy)  =  exp[ (-Um±iVm) y] 

This  can  be  accomplished  by  selection  of  the  proper  root  of 

and  /uTj\  If  we  illustrate  w^and  and  their  roots  in  the 
complex  plane,  we  have 


and  the  roots  of  Dn  +  iVm  arid  -Um  -iVn 

while  the  root.;  of  =  Um  -iVm  and  -Um  +iVm 
This  means  that  equations  (33)  and  (34)  become 

0  ym  =  2^  [Aiexp(i2irmx/A)+Biexp(“i2irrax/X)  ]exp[  (2tiio/A)  (-Umy-iVmy+if  At)  ] 

(35) 

B’ynt  -  S_i  [A_iexp(i27tmx/A)+B_]_exp(-i2TTmx/A)  ]exp[2Tim/A)  {-Umy+iVmy 

“if  At)  ]  (36) 


We  can  now  combine  the  particular  solutions  given  by  equations 
(35)  and  (36)  by  letting 


®1  =  A-1 
Ai  =  B_! 
Ei  =  E_i 


ym 


Adding  the  equations  together  gives: 

B'ym  =  Byxr«exP(_2irIQUmy/^)  (expi (im2it)  (ft+x/A-Vmy/X) 3  + 

expl(-im2ir)  (ft+x/A-Vmy/A)  ] } 
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or 


B*ym  =  Byiaexp(-2irittUniy/A) Cos(m2it/X)  (f Xt+x-Vmy)  (37) 

Before  going  any  further,  we  must  express  Um  and  Vra  in  terms  of 
y,  X,  f,  p’  and  m.  Since  Un  and  Vm  are  the  real  and  imaginary 
parts  of  /mi,  it  is  necessary  to  find  the  square  root  of  this 
complex  term, 

/ — ,  __  ..  iy2fX2.  x/2 
’'“l  =  d-f— 7--  -) 


=  /R{cos(8/2+irk)+isin(8/2+uk)  } 
k  =  0  and  1 

Since  we  previously  examined  all  the  roots  of  /ux*  and  /Si*,  we 
only  need  to  evaluate  this  case  for  k  =  0  to  determine  Um  and 
Vm- 


Thus, 


where 


l  =  »^{cosB/2+isin8/2} 


8  =  tan 


2fyX: 


Since 

cos 8/2  =  /l/zTl+coiir 


sin 8/2  =  /1/2  (l-cos8) 


The  equation  for  can  be  written  as 

4^  =  (1+Ml^iil) ^1+cosS’  +  i*  1-cos  6}’ 


p  •  -m 


It  can  be  noted  from  the  illustration  on  page  35  that 


ccsg  =  (I+li: 2y2^)  "i/2 
p  ’^m^ 
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We  can  now  form  a  series  of  solutions  which  will  satisfy  the 
boundary  condition  at  y  =  0  and  y  *  <*>.  By  letting 


00 


we  have  the  solution  for  the  By  component  of  the  flux  density 
in  the  pole  or 

By  «  BymexP("27rmUmy/AJcos^m2lI/x^  (f  Xt+x-Vmy)  (40) 

m=l,  2,  3 

In  this  case,  the  values  of  Vm  and  Vm  are  given  by  Equation  (38) 
and  (39)  and  Byia  equals  the  coefficients  for  the  cosine  terms 
in  a  series  which  give  the  "Y*  component  of  flux  density  at 
y  =  0.  We  can  also  form  a  series  of  solutions  for  the  B* 
component  of  flux  density  since 

V  -B  =  0  =  i§X.  + 

3x  3y 
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we  have 


f~  =  -|p^  =^^Bymexpilillexp(-2TrmUmy/X)[|unioos-2Y^(fXt  +  x-Vmy)  - 
m=l,2,3  Vmsin— ^(fXt  +  x-Vmy^ 


and  integrating  this  equation  with  respect  to  "x"  gives: 


Bx  =  y^ymgxp(-2yniUiny/X)jum  sin  ^  (fxt  +  x_Vjny) 

m=1,2»3  Vm  Cos  (fxt  +  x-Vmy )j 


Vm  Cos  (fxt  +  x 

A 

(y»t) 


(41  > 


The  question,  “Are  their  functions  such  as  (y,  t)  which  are 
a  part  of  the  complete  solution  for  the  B*  component  of  the 
flux  density?",  can  not  be  answered  since  we  have  no  boundary 
values  for  this  part  of  the  problem.  In  general,  the  solution 
given  by  Equation  (41)  for  (y,  t).  =  0  apparently  satisfies 
the  only  know  boundary  conditions  for  this  problem  which  are: 

(1)  Bjj  =  0  ,  y 

(2)  Bx  =  periodic  function  of  time 

(3)  V • B  =  0 

(4)  2-J&L  +  ilSi L  -  35x 

3x^  ay2  9t 


In  addition,  we  can  perform  two  integrations  which  shows  that 
the  flux  entering  the  pole  surface  is  conserved.  The  flux, 
0y,  entering  the  surface  when  ft  =  any  integer  is 


( x+1/r 

-k  Byl 


«tX 

1  dx 
f  y=Q 


f  x  +  ifi, 

-J  y  Byro 


This  flux  must  equal  Bx  wnen  By  =  0  or 


Cos  m27TXdx  =  XB-y^ 
X  2irm 


By  —  0 


,I^e*pC-2»mUmyA)dy=*|gSl 


The  lines  of  flux  in  the  pole  can  be  graphically  illustrated 
for  the  simple  case  when  m  =  1.  To  accomplish  this  task  we 
must  first  fc-rm  an  equation  for  the  slope  of  a  flux  line. 
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Since, 


dx 

ay 


lac 

Bv 


=  Ui  Sin  0  +  Un  Cos  0 
- Cos-0  - 


U-^tan  0  + 


where 

9  =  ~  (fXt  +  x-Viy)  =  f ( x , y , t ) 


and  since 


d© 


dx  +  i!  dy  +  dT 
9x  3y  3t 


we  will  let  ft  =  any  integer  and  dT  =  0.  Then 


dx 

ay  = 


d0  _  3f 
dy  ~  3y  ~ 
df 
3SE 


vi  *  vi 9  ♦  vi 

2 * 

X 


Thus 

dQ  -  2v  u-j  tan  9 
dy  x 

and  after  integrating 

2n  U-,  y  =  In  I  Sin  ol  +  constant 


Since  we  can  let  9o  =  0  at  the  point  where  y  =  0,  the  equa¬ 
tion  for  y  as  a  function  of  q  becomes 

2ttU-,  y  =  lr.  j Sin  0  | 

X  ISIiTffJ 


The  lines  of  flux  for  a  typical  simple  case  are  illustrated  in 
Figure  20.  The  only  parameter  which  one  needs  to  define  to 


room  temperature 


make  a  plot  line  the  one  shown  in  Figure  20  is  the  ratio  of 
vl/ul  since 

9  =  -  (Ml2>  (|jl) 


We  have  now  defined  the  rectangular  components  of  "B"  in  terms 
of  the  independent  variables  of  the  system.  The  next  task  is 
to  define  the  eddy-current  which  is  produced  by  the  flux. 

Our  solution  for  the  current  density  shall  be  based  upon  the 
following  considerations: 

(1)  Since  we  have  assumed  a  uniform  flux  in  the  "z" 
direction  which  does  not  vary  with  time,  we  shall 
also  assume  Jx  =  Jy  =  0.  Thus,  Equations  (16)  and 
(17)  are  of  no  importance  and  Equation  (18)  becomes: 

Jz  ,  52Jz  _  4 Tty  3 J z 
3x  3y2  ~  p*  3t 


and  Equations  (21),  (22),  and  (23)  reduce  to 

1_  3Bx  -  3 Jz 
p'  3t  3y 

1_  d  By  -  3 J  z 

P'  3t  "  3x 


(2)  Since  Bx  (x,  y,  t)  and  By  (x,  y,  t)  have  been  defined 
and  since  Jz  must  satisfy  the  same  partial  differen¬ 
tial  equation  as  By  and  Bjj,  we  can  obtain  Jz  from  the 
equations  listed  above.  However,  this  work  can  be 
avoided  by  using  Equation  (8).  For  this  case  we 
have 


,  T  3  By  3  Bx 
4TtyJ2  =  ^  -  3/ 


The  solution  for  Jz  is  obtained  by  forming  the  right-hand  side 
of  Equation  (42 ) 


4ityJz  = 


E[' 


<um-vm-l>  Sin  ^T2  (fxt  +  x~vmy>  * 


13=1,2,3 


4  U_V_  IB  Cos  llB 
m  m  X  X 


(fit  +  X-Vmy)J  exp(-2TrmUmy/A) 
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»«***  \'«*\  4  ‘W  <  Wfft* !- 


A 


2  2  f  llA  ^ 

Since  Un»-Vm-i  =  0  and  UmVm  =  -y- -  ,  this  equation  reduces  to 


fA 


JZ  = 


o'  /  i  Bym exp(-2irmUmy/A)Cos«~^fXt  +  X-Vmy) 
m=x^2 , 3 


J=  =  fr  % 


The  eddy-current  loss  per  unit  of  pole-face  is 


(43) 

(44) 


-HI, 


d+A 


p  * J^dxdy 


We  shall  evaluate  this  integral  for  the  case  when  ft 
integer.  Equation  (43)  becomes 


(45) 

any 


Jz  =  Sym  p  *  exp(-2TrmUJI^/A)cos51?iLXcosi2l!!lyjily  + 

m=l,2,3 


Bym  ^exp(-27rmUmy/x)sin?^^in2l2viny 


m=l ,  2  ,3 


If  we  form  the  product  of  p’  J2,  the  terms  in  the  series  would 
be  as  follows: 

B*  B  Lifiexpt - (mUm+nUn )2iry/A  ]Cos ZlIlMcos  L^Cos  ?.TOos  2,fmx  + 
J  3  P  A  A,  A  A 

BjayBny^-r— expC-  (mUm+nUr )  2iry/ A  }Sin  ilBZlliSin  glSgjB&in  linJfein  2nmx  + 

A  A  A  A 

2BmyBny^^exPf  ~ <roUm+nUn 5 2 »-¥/*  2  Co^lSgag Sirtl^SZ Co^irnxsin2inmt 
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where 


n  -  1 

and 

m  as  1, 

2, 

3 - k 

n  *  2 

and 

m  =  1 , 

2, 

3 - k 

n  =  3 

and 

m  =  1, 

2, 

3 - >k 

n  =  k 

and 

m  =  1, 

2, 

3 - k 

We  need  the  following  definite  integrals  to  evaluate  the 
integral  in  Equation  (45). 


m 


m 


n 


2 TfX  Cos 
A  " 


2ttx  Sin 
A 


2ttx  Sin 
A 


n  2ttx  dx 
A 


n  2irx  dx 
A 


m  2irx  dx 
A 


0  when  m  i  n 
1_  when  m  =  n 
2 


0  when  ra  i  n 
1  when  m  =  n 

T 


o 


With  these  definite  integrals,  all  integrals  with  sin  m  cos  n 

vanish  and  all  integrals  for  the  cases  when  m  is  not  equal 
to  n  vanish.  Hence,  integration  with  respect  to  ”x”  gives: 


P  = 


exp  ( -4irUmmy/ A )  dy 


(46) 


The  tasks  that  remain  are  (1)  to  determine  the  distribution 
of  flux  and  the  amplitude  of  flux  oscillation  on  the  surface 
of  the  pole,  (2)  to  determine  the  effective  resistivity  of 
the  material  when  the  pole  face  has  been  grooved,  (3)  compare 
the  experimental  results  with  results  obtained  from  theoretical 
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considerations,  and  (4)  to  determine  the  effective  permeability 
of  the  material. 


B.  amplitude  of  Flux  Oscillation  on  the  Surface  of  the  Poles 


A  solution  for  the  flux  distribution  on  the  surface  of  the 
pole  requires  a  solution  of  the  laplacian  differential  equation 


32F  ,  32F 
F‘x2+  r y2 


0 


An  analytic  solution  of  this  equation  can  be  obtained  by  con¬ 
formal  mapping  the  geometry  of  the  slots  and  the  opposing  rotor 
in  the  x-y-plane  onto  a  complex  plane.  One  then  constructs  a 
complex  potential  solution  in  order  to  determine  the  value  of 
field  vectors  in  the  complex  field.  Field  vectors  are  then 
expressed  in  terms  of  the  dimensions  of  the  x-y-plane.  In 
general,  this  approach  is  limited  to  idealized  configurations 
and  can  not  be  used  when  the  slots  are  closed  partially.  A 
numerical  or  graphical  solution  of  the  Laplacian  differential 
equation  is  required  when  the  configuration  of  the  slots  is 
complex . 

Since  the  magnetostatic  field  is  analogue  of  a  stationary 
temperature  field,  one  can  use  numerical  methods  developed  for 
heat  transfer  problems  to  determine  the  flux  distribution  on 
the  surface  of  the  pole.  In  this  case  the  correspondence  of 
quantities  in  the  fields  are: 


Quantity 

Potential  Function 


Basic  differential 
equation  of  the 
potential  function 

Characteristic  con’ 
stant  of  medium 

Potential  gradient 


Associated  field 
vector 

Flux  of  the  asso¬ 
ciated  vector 


Magnetostatic 

Field 


F^  magnetostatic 
potential 

=  0 


y  ,  absolute 
permeability 

H  =  -A  F,  mag¬ 
netizing  force 

B  =  uH  ,  Magnetic 
flux  density 

$  -  /  Bnds,  mag¬ 

netic  flux 


Stationary  Temp. 
Field 


T,  temperature 


A2T  =  0 


k,  thermal  conduc¬ 
tivity 

U  =  -a  T,  temperature 
gradient 

Q/A  =  kU  ,  heat 
power  flow  density 

Q  =  /  Q/A  dA,  heat 
power  flow 


% 
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The  problem  which  must  be  solved  is  to  determine  the 
magnetic  flux  density  distribution  on  the  surface  of  the  pole 
as  shown  below. 


fr 


The  equivalent  heat  transfer  problem  is  to  determine  the 
heat  power  flow  density  into  the  surface  of  pole  with  the 
surface  of  teeth  at  a  temperature  of  T,  and  with  the  surface 
of  the  pole  at  a  temperature  of  T2.  Tne  average  heat  power 
flow  density  on  the  surface  of  the  pole  would  be  subtracted 
from  the  value  obtained  at  each  point.  This  difference  would 
then  be  divided  by  the  average  heat  power  flow  density  to  give: 


Q/A-n/Aavg  _  B-Bavg 
Q/Aavg  “  feavg 


f (x/ x ) 


A  computer  program  is  available  to  IBM  7030  users  through 
the  Share  Library  which  will  solve  this  problem  by  a  numerical 
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approach.  The  program  is  called  TOSS  and  was  written  by  David 
Bagwell,  Union  Carbide  Nuclear  Company,  Oak  Ridge  Gaseous 
Diffusion  Plant,  Oak  Ridge,  Tennessee  (AEC  R  &  D  Report  No. 
K-1494).  TOSS  is  a  proven  digital  computer  program  for  finding 
the  transient  and/or  steady-state  temperature  distributions  of 
a  three-dimensional  irregular  body.  To  solve  the  problem  the 
body  is  divided  into  a  number  of  cells.  Each  cell  is  associated 
with  an  interior  node.  Surfaces  of  cells  which  lie  on  the 
boundary  of  the  body  are  associated  with  a  surface  node.  Con¬ 
ditions  external  to  the  system  which  affect  the  temperatures  of 
the  system  are  associated  with  boundary  nodes.  With  this 
nodal  description  of  a  system,  a  set  of  finite  difference 
equations  is  used  to  solve  heat  transfer  problems  that  are 
restrained  by 


k( 


32T 

3x'r 


+ 


32T)  x 
3Z2 


$C 


3T 


Internal  heat  generation  may  be  a  function  of  space  and 
time  and  boundary  temperature  may  be  a  function  of  time. 

A  nodal  network  which  described  the  space  between  the 
rotor  teeth  and  the  opposing  pole  face  in  test  fixture  was 
defined  for  an  air  gap  of  0.100  inches.  The  network  consisted 
of  302  interior  nodes,  50  surface  nodes,  and  2  boundary  nodes. 
The  first  boundary  node  represented  the  surface  temperature  of 
the  rotor  te^th  and  the  second  boundary  node  represented  the 
surface  temperature  of  the  pole.  A  steady-state  solution  of 
the  problem  was  obtained  which  gave  the  heat  power  flow  density 
at  18  points  along  one-half  of  the  tooth  pitch.  Figure  21  gives 
the  magnetic  flux  density  distribution  on  the  pole  surface 
that  was  obtained  with  TOSS. 

in  general,  the  magnetic  flux  density  distribution  on  the 
pole  surface  is  a  periodic  function  of  (x/x)  which  possesses 
certain  symmetric  properties  since 

f  (-x/  x)  =  f(*/x) 
and  j 


A  series  expansion  of  (Bavg-B)/Bavg  is  simple  since  it  is 
an  even  function  of  (x/X )  and  contains  no  sine  terms.  The 
coefficient  for  the  cosine  terms  can  be  obtained  by  making 
an  harmonic  analysis  of  the  periodic  function. 


Density  on  Surface  of  Poles  in  Experimental  Generator 


A  aeries  expansion  for  the  alternating  flux  density 
becomes 


3y  =  Bavg 

m=l , 2  }3 


Am  cos  m  x 


where  Am  equals  the  coefficients  obtained  from  the  harmonic 
analysis.  Letting  Bavg  Am  =  Bym#  this  equation  becomes 


By 


Bym  cos  m  22-  x 
m= 1 , 2  , 3  * 

Up  to  this  point  we  have  considered  "By"  on  the  surface  of  the 
pole  to  be  only  a  function  of  "x".  Since  "By"  varies  with 
time  at  any  point  on  the  surface  of  the  pole,  we  must  describe 
the  entire  function  in  order  to  establish  the  boundary  values 
for  the  problem. 


To  do  this  we  shall  consider  the  wave  shown  in  the  illus¬ 
tration  above  to  be  traveling  to  the  left  at  a  velocity  equal 
to  Xf.  We  shall  also  let  "x"  =  0  at  a  point  on  the  surface 
of  the  pole  where  "By"  has  its  maximum  value  and  the  product 
of  frequency  and  time  is  equal  to  any  integer.  Since  By  = 
f  (x,  t),  we  can  write 


dEL  =  dx  + dT 
3x  at 


and  since  dBy  =  0  when  -dx  =  XfdT  we  can  write 

Xf  =_£ 
ax  gt 

If  we  let  f ( x , t )  =  X*T 
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then 


t- 


f  „  X'  =  T*  =  •+  im2?tXf 
X  X  '  T 


and  the  solution  becomes 


3ym 


COS 


(f  At+X) 


(47) 


m=l,2,3 

Equation  (47)  gives  the  boundary  values  which  were  used  to 
obtain  By  (x,  y,  t)  3s  described  in  Section  IV-A. 


The  coefficients  for  the  cosine  terms  in  an  expansion  of 
the  flux  distribution  shown  in  Figure  21  are 


m 


1 

2 

3 

4 

5 


Arn  (coefficient ) 
0.0588 
-0.0087 
0.0008 


C . 0004 

0.0001 


Cosine  Term 
cos  2tfx/A 
cos  4?fX/X 
cos  6vx/\ 
cos  Stix/x 
cos  10tix/x 


The  coefficients  given  in  this  table  are  for  an  air  gap 
of  0.100  inches  or  for  the  air  gap  that  has  been  specified 
for  the  experimental  generator.  It  is  interesting  to  note 
that  for  this  particular  case,  over  98%  of  the  pole  face  loss 
in  the  generator  will  be  associated  with  the  fundamental  fre¬ 
quency  of  the  periodic  function.  The  loss  for  each  harmonic 
frequency  in  proportional  to 


Am  . 
mUm 


m 


n 

~  1  +  4f*ti* 

'V/2  +  i:oi 

__  |_  p“m 

"J  J 

TP. 


For  values  of  4f^u2A  *»  >>1 ,  the  loss  is  proportional  to  A2  .//in  • 

p  1  2m  m 

For  the  experimental  generator  and  for  a  0.100  inch  air  gap  in 


test  fixture. 

the  breakdown  of 

eddy-current  loss  is 

m 

/  /™ 

Proportion  of  Total  Loss 

1 

34.55  x  ltr4 

98.5% 

2 

53.5  x  10-6 

1.5% 

3 

37.0  x  10-8 

0 
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C.  Effective  Resistivity  of  a  Grooved  Pole  Face 


It  is  possible  to  reduce  the  eddy-current  loss  by  grooving 
the  pole  face  (reference  4).  The  grooves  must  be  placed  along 
the  face  of  the  pole  in  a  direction  parallel  to  the  motion  of 
the  rotor  so  they  cross  the  path  of  the  eddy-current s.  This 
increases  the  length  and  thus  the  resistance  of  the  eddy=*current 
path.  The  emf  from  end  to  end  of  the  pole  face  is  unaffected 
by  the  groove. 

The  resistance  of  the  path  under  the  grooves  is  directly 
proportional  to  the  depth  of  the  grooves  which  suggests  that 
the  deeper  the  grooves  the  greater  will  be  the  reduction  in 
loss.  However,  the  eddy-currents  can  take  an  alternate  path 
and  use  the  radial  face  of  each  groove  to  form  a  closed  loop. 

The  depth  of  the  groove  does  not  affect  the  resistance  of  the 
path  which  is  parallel  to  the  grooves. 


Since  I  =  1^  +  I2 

where  I2  =  current  flowing  under  the  grooves 

=  current  flowing  parallel  to  grooves 

The  total  resistance  offered  by  a  groove  is  ' 

Rg  =  Ri  +h2 

The  specific  resistance  of  the  path  under  the  groove 
1*2  =  {2d  +  w)p 

In  previous  work  it  has  been  shown  that  the  resistance  of  the 
path  parallel  to  the  grooves  is 

Ri  =  0.263  1  p 
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Thus 


Rg  +  (2d  +  w) »263x  p 
2d  *  w  +  . 2b3x 

The  resistance  Rg  is  in  series  with  the  natural  resistance  of 
the  pole  face  material.  Therefore, 

peff  .  *  (2d  +  w)»263xb  _  ,  (2d  +  w).263* 

p  "  1  (2d  +  w  +  .2b3AT  '  *  X  (2d  +  w  +  .263p 

g  (48) 

where 

d  =  depth  of  groove 
w  =  width  of  groove 
b  =  number  of  grooves 
1  =  length  of  pole  face 
X  =  tooth  pitch 
Xg  =  groove  pitch 

peff  =  effective  resistivity  of  pole  face  material 
p  -  natural  resistivity  of  material 

For  the  pole  face  loss  test  fixture, 

d  =  0.070  inches, 
w  =  0.007  inches, 
g  =  0.0325  inches, 

X  =  0.314  inches,  and 

Peff  =  2.63 
**T5  " 
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SECTION  V 


COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  RESULTS 
OBTAINED  FROM  THEORETICAL  CONSIDERATIONS 


The  eddy-current  loss  per  unit  area  of  the  pole  surface 
(watts  per  cm2)  versus  the  permeability  of  the  material  is 
given  in  Figure  22.  These  eddy-current  losses  were  calculated 
using  Equations  (38)  and  (46)  from  Section  IV-B  for  conditions 
which  correspond  to  a  racial  air  gap  of  0.100  inches  and  aver¬ 
age  flux  density  in  the  air  gap  of  56.5  kilolines  per  square 
inch  (8750  gauss).  Equation  (48)  was  used  to  calculate  the 
effective  resistivity  of  the  material  in  the  grooved  pole 
faces.  The  eddy-current  losses  given  by  Figure  22  are  those 
losses  associated  with  the  fundamental  frequency  of  the  periodic 
flux  variation  since  it  has  been  demonstrated  in  Section  IV-B 
that  the  losses  associated  with  the  harmonic  frequencies  are 
not  significant  for  this  particular  case. 

Figure  22  also  shows  the  pole  face  losses  that  were  derived 
from  the  experimental  data  obtained  in  the  pole  face  loss  test 
fixture.  A  comparison  of  losses  derived  from  the  experimental 
data  and  the  losses  based  upon  theoretical  consideration  indi¬ 
cates  that  the  effective  permeability  of  the  material  is  extreme¬ 
ly  low.  In  fact,  the  experimental  loss  at  8000  cps  is  greater 
than  the  losses  calculated  for  a  material  which  has  3  perme¬ 
ability  of  one.  Since  the  effective  permeability  appears  to 
be  low,  ve  must  determine  the  maximum  induction  in  the  pole  to 
determine  if  the  material  is  saturated.  The  maximum  flux 
density  occurs  on  the  surface  of  the  pole  or  when  y  =  C. 
Therefore, 

Bmax  =  [(Bavg  +  By)  2  *v  £x2]-5 

=  BaVg  [1  +  A^U^  +  2AX  (Cose  +  .5AiUiVxsin26)  3  *5 


where 

9  =~  (f,\t+x) 

Figure  23  give  B^x  as  a  function  of  Q  for  a  fundamental 
frequency  of  24,000  cps  and  for  a  value  of  U^Vi  and  A^  which 
corresponds  with  the  loss  observed  in  the  test  fixture.  The 
maximum  flux  density  for  this  case  was  only  6.2%  greater  than 
the  average  flux  density  in  the  air  gap.  Figure  24  shows  the 
maximum  and  minumum  flux  density  at  the  surface  for  p  (effective) 
=  10.5  on  a  normal  magnetization  curve  for  the  H-ll  material. 
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DIMENSIONLESS  DISTANCE  ALONG  FACE  OF  POLE  (ELECTRICAL  DEGREES) 
Figure  23.  Maximum  Flux  Density  in  Pole  Faces  at  24,000  cps 
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H,  OERSTEDS 


Figure  24.  Maximum  and  Minimum  Flux  Density  in  Pole  Face  at 
24,000  cps  Graphically  Illustrated  on  a  Normal 
Magnetization  Curve  for  H-1I  Steel 
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It  can  be  noted  that  the  maximum  induction  is  well  below  the 
point  of  saturation.  Thus,  the  low  value  of  p  (effective)  can 
not  be  associated  with  points  of  saturation  within  the  pole  face. 
The  low  value  of  p  (effective)  may  also  be  due  to  the  fact  that 
a  significant  portion  of  the  tooth  ripple  loss  is  produced  by 
hysteresis  loss.  Minor  hysteresis  loops  are  produced  by  the 
varying  field  which  are  superposed  upon  the  normal  magnetization 
curve  as  shown  in  Figure  24.  A  different  loop  would  exist  for 
each  point  in  the  pole.  The  loop  would  represent  the  variat-on 
of  B  and  H  with  time  for  a  constant  value  of  X  and  y.  The 
hysteresis  loss  in  watts  per  cubic  cm  at  a  point  is 


p 

=  f  xlO-7 

f  B(t ) 

X=C1  8  *  ' 

Jc 

y=c2 

d[Hi't>] 

x=c^ 

y=c2 


and  the  total  hysteresis  loss  per  unit  of  surface  area  would  be 


fxlO 
8  is 


-7  ( 


JU 


B(t) 


X=Cj^ 

y=c 


d[H(t)]j  dy 
Jx=c 
y=0 


to 00 


In  general,  the  area  within  a  hysteresis  loop  can  only  be 
evaluated  experimentally.  However,  some  conclusions  can  be 
obtained  if  one  assumes  the  following: 


(a)  the  area  within  the  hysteresis  loop  at  y  =  0  and 
X  =  0  shall  be  assumed  to  be  equal  to 

/  Bdh  =  2A^  Bavg  Havg 
c 

wh.ich  is  probably  several  times  greater  than  the 
actual  case. 

(b)  the  area  within  the  hysteresis  loops  at  X  =  0  shall 
be  assumed  to  decrease  with  y  in  accordance  with  the 
relation 

Area  =  (Area)y_0exp(-21rU1y/X) 
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This  relationship  is  based  upon  the  assumption  that  the 
area  within  a  hysteresis  loop  is  proportional  to  the  amplitude 
of  cyclically  changing  flux  density.  In  the  actual  case  the 
area  would  also  be  proportional  to  the  amplitude  of  the 
cyclically  changing  field- strength.  Thus, 

Area  =  (Area)y_Qexp(-4iTUiy/X) 

is  probably  more  correct  than  the  relationship  assumed  for  the 
calculation  presented  in  this  report. 

'.Jith  these  assumptions  hysteresis  loss  per  unit  of  surface 
area  becomes 


P  u  =  A-]  BdVj 


Havg  X  fxlO* 
8ir2U-, 


(watt /cm2) 


Equation  x49)  gives  an  extravagant  value  for  the  hysteresis 
loss.  Thus,  if  the  value  of  hysteresis  loss  obtained  from  this 
equation  is  small  compared  with  the  calculated  eddy-current  loss, 
the  inclusion  of  hysteresis  loss  would  not  be  essential.  Figure 
25  gives  pole  face  losses  with  and  without  the  inclusion  of 
hysteresis  loss.  In  general,  this  comparison  shows  tha»_ 
hysteresis  loss  probably  is  not  a  significant  part  of  the  total 
tooth  ripple  loss  for  these  operating  conditions.  The  inclusion 
of  hysteresis  loss  means  that  the  effective  permeability  for 
eddy-current  calculations  would  have  to  be  increased  slightly  to 
provide  for  a  small  decrease  of  eddy-current  loss  that  is 
required  for  the  same  combined  losses.  Figure  25  shows  the 
eddy-current  loss  is  predominant  over  hysteresis  3oss  at  high 
tooth  ripple  frequencies.  Therefore,  the  inclusion  of  hystere¬ 
sis  loss  is  not  essential  for  the  operating  conditions  which 
establish  the  maximum  heat  load  for  the  rotor  cooling  system. 


The  only  reasonable  explanations  for  the  low  values  of 
effective  permeability  that  were  derived  from  the  experimental 
data  are:  (1)  that  H-ll  steel  had  a  low  reversible  perme¬ 
ability,  and  (2)  that  he  incremental  permeabilities  associated 
with  the  minor  hysteresis  loops  were  also  low.  Figure  26 
illustrates  how  minor  hysteresis  loops  under  various  conditions 
in  iron  appear  on  a  normal  magnetization  curve.  The  reversible 
permeability,  yr  ,  for  value  of  Ej,  and  B]-  lying  on  a  normal 
magnetization  curve  can  be  determined  by  two  ways  (reference  5  ) . 
The  first  procedure  is  illustrated  on  Figure  26b.  In  this  case 
the  biasing  field  strength  is  held  constant  and  the  incremental 
permeability, w A  =AB/£H,  is  determined  for  an  increasing  AH. 

The  reversible  is  then  obtained  by  extrapolating  to  AH  =  0. 
The  second  procedure  is  illustrated  on  Figure  26c.  For  this 
case,  a  steady  field  is  applied  in  one  direction,  and  the 
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Fxgure  25.  Pole  Face  Losses  with  the  Inclusion  of  Hysteresis 

Loss 
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alternating  field  is  applied  first  in  the  opposite  direction. 
The  incremental  permeability,  y^  =  AB/AH,  is  determined  for  an 
increasing  AH  and  the  reversible  permeability  is  determined  by 
extrapolating  to  AH  =  0. 

The  reversible  permeability  of  the  pole  face  material  is 
an  important  parameter  since  it  is  the  applicable  value  of  per- 
meavility  to  be  used  in  the  pole  face  loss  equations  when  the 
value  of  AH  or  Alt  approach  zero.  In  the  actual  case  AH  and  Ab 
approach  zero  at  points  under  the  surface.  Therefore,  the  re¬ 
versible  permeability  would  be  the  exact  value  of  permeability 
to  be  used  for  eddy-current  losses  at  these  points.  At  a  point 
in  the  pole  where  AB  has  a  finite  value,  the  permeability  ap¬ 
proaches  the  incremental  permeability  of  the  minor  hysteresis 
loop  which  is  produced  at  this  point.  The  incremental  perme¬ 
ability  represents  the  average  change  of  B  with  respect  to 
changes  of  H  along  the  hysteresis  loop. 


The  applicable  value  for  the  permeability  for  calculating 
pole  face  losses  can  be  determined  by  a  numerical  integration. 
Equation  46  gives  the  pole  face  loss  as 


P  =  (X2f?2p  •) 


ST  B2ym  exP 

0  m=l ,2,3 


(-4Trl^ny/A)dy 


The  integrands  of  this  equation  are  equal  to  the  square  of  the 
y  component  of  flux  density  for  each  harmonic  frequency  along 
a  path  where  /ft  4-  x/A-Vmy/A)  equals  any  integer.  If  we  con¬ 
sider  each  harmonic  frequency  separately,  the  loss  for  each 
frequency  becomes 


Pm 


(a zf 2/2p  * ) 


2 

By  dy 


Since 


By  =  Bym  exp  {— 2TimtImy/A) 
dy  =  (A/2nifegnBy)d  By 

we  can  express  the  loss  for  each  frequency  as 
Pm  =  ( A3f  3/8itmp  1 }  [Bym 

J  V m 

o 

In  order  to  simplify  the  numerical  integration  of  this  equation, 
we  shall  assume  that  4f2u2A*I/p  ,2m2>>l  and  it  can  be  noted  from 
Equation  38  that  Um  becomes  equal  to  A/fy/p 'm.  Thus, 


1 
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(50) 


Pm  = 


(X2f3/V3  v/PE) 


The  integrand,  l//y",  in  the  equation  shall  be  considered  to  be  a 
variable  quantity  and  equal  to  CH//y7.  By  is  equal  to  the 
amplitude  of  the  y  component  of  the  flux  density  for  a  constant 
value  of  "y" .  If  we  consider  a  constant  value  of  "y**,  the  flux 
density  has  a  sinusoidal  variation  with  "x"  or  time  which 
produces  minor  hysteresis  loops  as  illustrated  by  Figure  26C. 

The  AB  associated  with  each  loop  is  equal  to  2By.  The  biasing 
induction,  will  be  constant  and  equal  to  the  average  flux 
density  in  tne  pole.  A  different  loop  exists  for  each  value  of 
nyn .  Since  ma  varies  with  the  amplitude  of  the  loop,  a  numerical 
integration  of  Equation  (50)  is  required  as  illustrated  in 
Figure  27.  The  parameter  "Cfl"  is  a  correction  for  the  fact 
that  the  slope  at  a  point  on  a  minor  hysteresis  loop  is  not 
equal  to  y& .  The  previous  derivations  were  based  upon  the 
assumption  that  y  =  yA  and  that  the  loss  per  unit  volume  was 
proportional  to  ByVy^ .  A  largo  portion  of  a  minor  hysteresis 
loop  has  a  slope  which  is  less  than  at  points  where  Byz  is 
large.  Therefore,  a  significant  error  was  introduced  when 
(AB/2)  2  (i//yX)  cos1 2 3 4 52t?x/A  was  integrated  with  respect  to  "x".  The 
factor  required  to  correct  this  error  is 


CH  = 


/•2TT 

/  {J 
J  o 


J ^  (1  //y")  cos2  {2irx/A)d{27rx/A) 
-Q&:  I  cos2  (2irx/A)d(2irx/.A) 

,vtJ, 


where  y  =  the  actual  slope  at  a  point  on  the  hysteresis  loop 
and  B&  =  AB/2.  The  evaluation  of  the  integral  requires  a  stand¬ 
ardization  of  the  hysteresis  loops.  The  specific  assumptions 
made  regarding  the  shape  of  the  loop  as  shown  on  page  64  are: 


(1)  the  slope  of  the  hysteresis  lcop  is  a  linear  function 
of  the  field  strength 

(2)  the  slope  of  the  hysteresis  loop  for  decreasing  in¬ 
duction  at  By/BA  s  1  is  equal  to  the  reversible 
permeability 

(3)  the  slope  of  the  hysteresis  loop  for  increasing  in¬ 
duction  at  By/BA  =  -1  is  equal  to  the  reversible 
permeability 

(4)  Cos“l{By/BA)  =  2’rx/A 

(5)  for  By/B^  ~  1  H/Ha  =  1  and  for  By/B a  «  -1  H/Ha  =  -1 
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Figure  27.  Effective  Permeability  of  11-11  Steel  Based  Upon 

Theoretical  Considerations 


63 


with  these  assumptions 
al  =  a2  =  pa 

®1  “  c2  - Jur-PAJHA/2 

and  for  the  upper  portion  of  the  loop, 

0  <  x/A  <1/2  |  P/UA  *  1  +  yA)H/yAHA  (5i) 

/By/BA  =  ((yr-  yA)/2UA)  (H/Ha}2  +  H/HA-(yr-JlA)/2uA 

(52) 

For  the  DOttcm  portion  of  the  loop  the  equations  become 

1/2  1  x/A  <  1  |  n/yA  *  1  ~  {ur-HA>H/yAHA  (53) 

I  By/BA  =  -  ((vz-vi)/2Va)  (H/Ha)2  +  H/Ha  +  (vr-uA)/2pA 

(54) 
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Since  the  correction  factor  for  the  upper  portion  of  the  loop 
will  be  equal  to  the  correction  factor  for  the  lower  portion, 
we  only  need  to  determine  the  factor  for  one-half  of  the  loop. 
The  task  that  remains  is  to  form  the  integrals  required  to 
evaluate  the  correction  factor.  The  correction  factor  can  be 
expressed  as 


CH  =  ± 


^  J  (By/BA)  2d(27rx/X) 

0 


Since  and  p^are  not  a  function  of  ”x"  and 


By/BA  =  cos  2ux/A 
By  letting 


CH  ==  2/ 


'  / 
J  \ 


-1 


f  {By/BA)d(By/BA) 


-2/'l 


/i^77(By/BA)  ^d(2Trx/A) 


we  can  form  the  integral  on  the  left, 
because 

cos-1 (By/BA)  =  2n  x/A 

-  d(2™A> 

and  from  equations  51  and  52  it  can  be  shown  that 

Therefore,  the  integral  for  the  correction  becomes 


CH  =  2/ 


/“  -(By/BA)2d(By/BA) . 

j  /l-(By/B A)  2  |  1  +  ^Pr~P&^2t2^Mr~PA^  Byj.1 '  ** 


if  we  let  By/BA  =  Cose  and  (pA-ur)/uA  =s, 
Cg  ~  2/n  j  cos2ede 

J  o 


{l+e2-2eoos6 } 1 ' ^ 


(55) 
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The  correction  factor  for  various  values  of  e  was  determined 
by  Gaussian  quadrature  integration  of  Equation  (55)  as  shown 
by  Figure  28.  The  applicable  value  of  permeability  for  pole 
face  loss  calculations  can  now  be  determined  with  the  correc¬ 
tions  given  by  Figure  28  and  with  the  numerical  integration 
shown  on  Figure  27 .  The  material  properties  required  for  this 
integration  are: 

. (1)  the  reversible  permeability  of  the  material  at  a 

steady  induction  equal  to  the  average  flux  density 
in  the  poles  and, 

(2)  the  incremental  permeabilities  of  the  minor  hysteresis 
loops  which  occur  in  the  poles. 

Since  these  quantities  are  intrinsic  and  measurable  properties 
of  the  material  the  effective  permeability  is  a  predictable 
quantity  once  the  flux  distribution  on  the  surface  of  the  pole 
is  known .  Unfortunately,  the  reversible  and  incremental  per¬ 
meabilities  of  H-ll  steel  have  not  been  measured  and  it  is  there¬ 
fore  necessary  to  derive  a  value  for  the  effective  permeability 
of  H-ll  steel  from  theoretical  considerations.  When  the  ampli¬ 
tude  of  the  alternating  field,  HA,  is  increased  the  incremental 
permeability  is  also  generally  increased.  For  small  amplitudes 
this  increase  is  linear  and  by  analogy  with  Rayleigh’s  law, 
we  can  assume  that 

ua  =  Pr  + 


and 


ba  =  WrH&  +  vaha2 

When  the  magnetization  is  sufficiently  low  Rayleigh's  law  may 
be  used  for  the  initial  portion  of  the  normal  magnetization 
curve  or 

v  =  p0  +  VH 


and 


3  =  p0H  +  VH2 

For  H-ll  steel  the  initial  portion  of  the  normal  magnetization 
curve  gives 

y  o  -  20 
V  =  6 
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Ebinger,  reference  5,  has  measured  V&  for  a  variety  jii  materials 
and  plotted  V&/V  vs.  (B-H)/Bs*  He  found  that  the  points  lie 
reasonably  close  to  a  single  curve.  For  a  biasing  induction  of 
8.75  kilogausses  on  H-ll  steel  the  ratio  of  Va/V  from  Ebinger *s 
curve  is  approximately  0.175  which  means  that  value  of  Va  for 
H-ll  steel  at  this  induction  is  approximately  one.  Brown, 
reference  5,  has  shown  for  material  with  no  direction  of  easy 
magnetization  and  with  anisotropic  domains  that  the  ratio  of 
ur/vo  equals  approximately  0.15  for  a  proportionate  induction, 
(B-H)/Bs,  of  42%  which  means  that  Pr  for  H-ll  steel  at  an  in¬ 
duction  of  8.75  kilogausses  would  be  approximately  equal  to 
three. 

Thus,  the  equation  for  the  incremental  permeability  becomes 
HA  =  3  +  1.05  Ha 

and 

Ba  =  3Ha  +  1.05H| 

We  can  now  express  pa  *  vr/vA*  Ba ,  and  other  parameters  of 
interest  in  terms  of  Ha  as  shown  below. 


H 

oersteds 

0 

“a 

3.0 

ba 

gausses 

0 

vx/v& 

1.0 

CH 

1.0 

CH//^ 

0.577 

Bi 

(gausses) 

0 

2 

2 

5,1 

10.2 

0.588 

1.038 

0.460 

0.000104 

x 

106 

4 

7.2 

28.8 

0.416 

1.081 

0.402 

0.000830 

X 

106 

6 

9.3 

55.8 

0.322 

1.115 

0.368 

0.00311 

X 

106 

8 

11.4 

91.2 

0.263 

1.140 

0.331 

0.00781 

X 

106 

10 

13.5 

135.0 

0.222 

1.160 

0.316 

0.0182 

X 

1C6 

12 

15.6 

187.0 

0.1922 

1.180 

0.299 

0.035 

X 

lOf 

14 

17.7 

248.0 

0.1694 

1.195 

0.283 

0.0615 

X 

106 

16 

19.8 

317.0 

0.1513 

1.206 

0.271 

0.1003 

X 

106 

18 

21.9 

394.0 

0.137 

1.216 

0.260 

0.1552 

X 

lof 

20 

24.0 

480.0 

0.125 

1.226 

0.250 

0.230 

X 

106 

22 

Since 

26.1 

574.0 

0.115 

1.234 

0.242 

0.329 

X 

10® 

i  r By®  c 

/jT  effective=  ^/Bym  /  /=  ^^BA^ 

o 

the  applicable  value  of  permeability  can  be  determined  by 
numerical  integration  of  the  integral  as  shown  in  Figure  27. 
This  integration  gives 
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l//jT  effective  =  0.274 


* 


for  Bym  -  515  gausses.  Thus, 


jj  effective  =  14.3 

Although  the  above  analysis  may  not  be  strictly  accurate 
from,  a  numerical  standpoint  due  to  the  admission  ~>f  approxima¬ 
tions,  it  shows  that  the  effective  permeability  of  H-lI  steel 
is  low. 

It  also  shows  that  the  low  permeability  value  required  to 
obtain  agreement  between  experimental  and  analytical  results  is 
reasonable . 

The  eddy-current  losses  given  by  Figures  22  and  25  are 
based  upon  conditions  pertinent  to  the  experimental  generator 
except  for  resistivity  of  the  material.  The  losses  given  by 
Figures  22  and  25  were  calculated  using  resistivity  which  corre¬ 
sponds  to  room  temperature  conditions.  Figure  29  gives  the 
eddy-current  loss  with  a  correction  for  a  change  in  resistivity 
based  upon  the  estimated  operating  temperature  of  the  pole  face 
in  the  generator. 
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SECTION  VI 


POLE  FACE  LOSSES  CAUSED  BY  ARMATURE  REACTION 


The  eddy-currents  produced  by  the  stray  rotating  fields 
of  the  armature  reaction  were  calculated  by  using  the  equations 
from  Reference  2.  These  equations  were  published  in  1955  by 
G.  Barello  and  his  derivations  shall  be  taken  for  granted.  His 
results  give: 


£ 


JL.T  r 1  t  fei2  3/2  3/2 

5.65  lQQj  LlQOJ  L  n  J  (nil)  fi 

'\ivof  sinhimS  +  T-;m  cosh  irn6  2  +  T^m  coshimS 
l  T^l  yrivn 


where : 


W  =  Loss  in  watts  per  square  meter  of  pole  surface 
N  =  Number  of  conductors  in  series  per  pole  per  phase 
I  =  Current  per  phase  in  amperes 

fjL=  Nominal  frequency  of  currents  in  windings,  the  base 
frequency  of  the  machine  in  cps 

n  =  The  order  of  the  harmonic  ~»eing  considered 

fen  =  Kp  Kd  for  the  nth  harmonic  being  considered 

pr  =  The  relative  permeability  of  the  rotor  material 

p  =  Resistivity  in  ohm-meters 

Tj  =s  Pole  pitch  of  the  machine  in  meters 


6  =  Air  gap  in  meters 
m  =  /iTvi(ntl  )f  ]_/p 


p  u  oPr 

u0=  1.256  x  1CT6 


in  meter3 
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Barello  found  that  a  value  of  yr  =  2000  was  a  suitable 
value  for  the  relative  permeability  of  the  rotor  material.  In 
general ,  this  value  gave  good  agreement  between  calculated 
losses  and  the  losses  measured  on  tests  with  similar  machines. 

The  calculated  eddy-current  losses  produced  by  the  stray 
rotating  fields  of  the  armature  reaction  for  the  experimental 
generator  are  presented  in  Table  I.  For  this  case  the  losses 
are  20.9  kilowatt  per  square  meter.  Since  the  total  area 
of  the  pole  faces  in  the  generator  is  0.050  square  meters, 
the  losses  due  to  the  stray  rotating  field  of  the  armature 
reaction  will  be  1050  watts. 

Unfortunately,  the  original  calculated  value  for  the 
losses  due  to  armature  reaction  was  210  watts  and  the  rotor 
cooling  system  was  designed  for  a  combined  value  of  pole  face 
losses  equal  to  approximately  2000  wats..  Thus,  with  1050 
watts  of  loss  due  to  armature  reaction  and  1793  watts  loss 
due  to  tuffs  of  flux  from  the  stator  teeth,  over-heating  of 
the  rotor  will  occur  if  the  generator  is  operated  at  design 
conditions.  This  conclusion  is  supported  by  the  fact  that 
the  calculated  loss  due  to  armature  reaction  was  based  upon 
a  relative  permeability,  ur,  of  2000  for  the  rotor  material. 
This  value  had  an  empirical  derivation  based  upon  test  data 
from  average  machines  and  it  is  reasonable  to  assume  the 
applicable  value  was  influenced  by  the  actual  permeability 
of  the  rotor  material.  Since  the  material  in  the  rotor  of 
the  experimental  generator  has  an  unusually  low  permeability, 
it  is  problematical  that  the  losses  due  to  armature  reaction 
will  even  be  higher  than  the  calculated  value  of  1050  watts 
based  upon  Barello' s  equation.  In  any  case,  the  loss  due 
to  armature  reaction  constitute  a  large  portion  of  the  total 
estimated  pole  face  losses. 


TABLE  I.  Pole  Face  Losses  Caused  by  Stray  Rotating  Fields  of 
Armature  Reaction  Losses  by  Bareilo  Method 
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SECTION  VII 


SUMMARY  AND  CONCLUSIONS 


The  pole  face  losses  caused  by  the  alternating  flux 
density  that  is  produced  by  the  slots  and  teeth  of  the  arma¬ 
ture  for  the  467  leva  experimental  generator  were  experimentally 
evaluated  by  using  an  eddy-current  dyanamometer .  Design  factors 
which  influence  pole  face  losses  in  the  generator  were  duplicated 
in  the  test  fixture  as  close  as  possible.  Centrifugal  stress 
limited  the  rotor  speed  in  the  test  fixture  to  one-half  the 
speed  of  the  generator.  It  was,  therefore,  necessary  to  extra¬ 
polate  the  data  from  a  tooth  ripple  frequency  of  24,000  cps  to 
48,000  cps  to  determine  the  face  losses  in  the  generator.  No 
attempt  was  made  to  control  the  temperature  of  the  pole  faces. 
Thus,  the  test  results  had  to  be  corrected  for  a  change  of 
resistivity  since  the  pole  faces  will  operate  at  a  much  higher 
temperature  in  the  generator.  The  pole  face  losses  due  to  tooth 
flux  ripple  for  the  generator  as  derived  from  the  experimental 
data  are  1783  +  590  watts  for  pole  face  temperature  of  950°F. 

An  analysis  of  the  losses  associated  with  tooth  flux  ripple 
based  upon  theoretical  considerations  was  made.  A  procedure 
for  dealing  quantitatively  with  the  magnetic  flux  density 
distribution  at  the  surface  of  the  pole  face  was  outlined.  A 
knowledge  of  the  distribution  of  flux  density  on  the  surface 
is  essential,  because  it  is  the  non-uniformity  of  this  flux 
which  causes  all  tooth-ripple  losses.  Harmonic  frequencies  and 
losses  due  to  each  harmonic  frequency  were  considered.  This 
analysis  showed  that  98 . 5%  of  the  tooth-ripple  losses  in  the 
experimental  generator  are  associated  with  the  fundamental 
frequency  of  the  cyclically  changing  induction  in  the  pole 
faces.  The  amplitude  of  the  cyclically  changing  flux  density 
for  the  fundamental  frequency  was  found  to  be  5.88%  of  the 
average  flux  density  in  the  radial  gap.  Solid  and  grooved  pole 
faces  were  considered  in  the  analysis.  The  losses  in  grooved 
pole  faces  were  calculated  by  using  an  effective  resistivity  for 
the  pole  face  material  which  was  based  upon  the  increase  in  the 
resistance  of  the  eddy-current  path  caused  by  the  grooves. 

The  resistivity  of  the  pole  face  material  was  effectively 
increased  by  2. 63 -folds.  Thus,  grooving  the  pole  faces  re¬ 
duced  the  pole  face  losses  by  approximately  38.5%.  Hysteresis 
losses  in  the  pole  faces  could  not  be  separated  from  the  other 
losses  since  no  theoretical  relationship  exists  between  the 
amplitude  of  flux  density  oscillation  and  the  area  within  a 
hysteresis  loop.  In  order  to  determine  if  hysteresis  losses 
were  of  importance  a  relationship  between  the  area  of  a  hysters- 
sis  loop  and  the  amplitude  of  flux  density  oscillation  was 
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assumed  which  would  give  an  extravagant  value  for  the  hysteresis 
loss.  Even  for  these  conditions  it  was  found  that  hysteresis 
loss  was  of  minor  importance  and  that  eddy-current  loss  was 
predominant  over  hysteresis  loss. 

The  theoretical  equations  for  the  calculations  of  the 
eddy-current  losses  contain  a  parameter  which  is  equal  to  the 
permeability  of  the  pole  face  material .  The  derivations  of 
these  equations  were  based  upon  the  assumption  that  the  perme¬ 
ability  of  the  material  is  a  constant.  In  the  actual  case  the 
permeability  of  the  material  is  a  variable  and  an  applicable 
value  of  this  parameter  must  be  used  to  obtain  a  reasonable 
estimation  of  the  eddy-current  losses  from  theoretical  consider¬ 
ation.  The  experimental  results  present  in  this  report  demon¬ 
strate  that  the  applicable  value  of  permeability  is  associated 
with  the  incremental  permeability  of  the  minor  hysteresis  loops 
that  are  produced  in  the  pole  face  material.  In  general,  the 
incremental  permeability  of  a  minor  hysteresis  loop  in  a  parti¬ 
cular  material  can  vary  widely  since  the  slope  of  a  loop  usually 
increases  when  the  amplitude  of  the  oscillating  induction  increas 
es.  A  procedure  for  dealing  with  a  varying  incremental  perme¬ 
ability  was  outlined.  The  applicable  value  can  be  obtained  by 
a  numerical  integration  which  also  considers  the  varying  perme¬ 
ability  around  each  minor  hysteresis  loop. 

The  eddy-current  loss  produced  by  the  stray  rotating  fields 
of  the  armature  reaction  were  calculated  using  Barello's  equa¬ 
tions.  It  was  concluded  from  these  calculations  that  the  rotat¬ 
ing  fields  of  armature  reaction  will  produce  a  loss  of  1050  watts 
Unfortunately,  the  rotor  cooling  system  was  designed  for  only 
210  watts  of  loss  due  to  armature  reaction.  This  dilemma  results 
from  inadvertant  error  in  the  original  calculations.  It  is 
therefore  generally  concluded  that  overheating  of  rotor  would 
occur  if  the  generator  was  operated  at  full  load  with  a  575°F 
coolant  inlet  temperature.  This  particular  dilemma  was  not 
resolved  by  redesigning  the  cooling  system  since  the  redirected 
test  plans  for  the  generator  allows  a  lovr  temperature  coolant 
to  be  used.  Nevertheless,  it  is  a  logical  conclusion  that 
the  eddy-current  loss  produced  by  the  rotating  fields  of  the 
armature  reaction  will  constitute  a  large  port.^n  of  the  rotor 
face  losses  in  this  type  of  generator. 

The  design  of  an  electromagnetic  generator  with  a  high-speed 
high -temperature  solid  rotor  must  cope  with  the  problem  of  re¬ 
ducing  the  pole  face  losses  if  temperature  is  to  be  attained 
which  do  not  impose  serious  design  limitations  upon  rotor 
materials.  There  are  certain  parameters  which  can  be  varied 
to  reduce  the  rotor  pole  face  losses.  Some  of  the  parameters 
which  have  a  significant  influence  upon  the  pole  face  losses 
are  the  radial  gap,  the  ratio  of  the  width  of  the  slot  opening 
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and  tooth  pitch, tooth  pitch,  pole  pitch,  and  the  direct  axis 
synchronous  reactance  of  the  machine.  In  general,  these  para¬ 
meters  can  be  varied  to  reduce  the  amplitude  of  flux  density 
variations  on  the  surface  of  the  pole  faces.  However,  serious 
performance  and  weight  penalties  can  be  incurred  in  the  overall 
design  when  these  parameters  are  optimized  only  from  the  stand¬ 
point  to  pole  face  losses.  The  designer  is  therefore  obliged 
to  look  for  other  methods  for  reducing  pole  face  losses.  The 
most  promising  solutions  of  the  dilemma  ’would  be  to  combine  a 
material  which  offers  excellent  properties  for  the  reduction  of 
eddy-current  losses  in  the  pole  faces  with  a  material  that  offers 
excellent  creep  strength  and  d-c  magnetic  properties  for  a 
solid  rotor.  Pole  face  losses  are  reduced  by  selection  of  a 
material  for  the  pole  tips  which  has  a  high  resistivity  and  a 
high  reversible  permeability.  Unfortunately,  materials  which 
have  the  best  combination  of  d-c  magnetic  and  strength  pro¬ 
perties  for  a  high  speed  and  high  temperature  rotor  have 
inferior  properties  from  the  standpoint  of  pole  face  losses 
For  example,  the  pole  face  losses  in  Westinghouse  NIVCQ  alloys 
will  be  approximately  2  to  3 -folds  greater  than  the  pole  face 
losses  in  H-ll  steel.  The  losses  in  H-ll  steel  are  approximately 
4-folds  greater  than  the  losses  in  soft  iron.  The  increase  in 
heating  in  the  pole  faces  may  more  than  offset  the  advantage 
offered  by  the  additional  high-temperature  creep  strength  of 
a  particular  material. 

Generally  a  material  which  has  a  high  normal  permeability 
at  a  particular  value  of  induction  also  has  high  incremental 
and  reversible  permeabilities  at  this  level  of  induction.  The 
normal  permeability  of  the  cobalt-iron  alloys  at  inductions  up 
to  10  kilogausses  are  relatively  high.  Thus,  these  materials 
and  probably  others  are  candidates  for  pole  face  materials.  It 
should  be  noted  that  a  material  of  much  lower  creep  strength 
can  be  used  on  the  pole  tips  since  the  stress  is  much  lower  on 
the  outer  radius  of  the  rotor  than  at  the  root  of  the  pole 
bodies  or  in  the  rotor  core.  In  addition  the  creep  strain  may 
be  allowed  to  approach  1  to  1-1/2%  since  the  length  of  this 
material  will  be  small .  A  creep  strain  over  1/2%  for  Other 
parts  of  the  rotor  may  be  considered  excessive  since  this 
amount  of  strain  would  cause  a  large  change  of  the  radial  gap 
in  the  generator . 

A  large  increase  in  resistance  tc  the  flow  of  eddy-current 
can  be  obtained  with  laminated  pole  tips.  Consequently ,  arrange¬ 
ment  for  attaching  the  pole  tips  should  permit  a  laminated  con¬ 
struction  to  be  used.  Eddy-current  losses  can  be  decreased  by 
approximately  2-foius  by  using  laminations  in  place  of  grooved 
pole  faces. 
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The  overall  decrease  of  eddy-current  losses  with  a  laminated 
material  which  has  been  chosen  for  its  high  a-c  magnetic  perme¬ 
ability  can  approach  16-folds  when  compared  with  the  losses 
associated  with  H-ll  steel  or  other  materials  that  have  been 
selected  for  their  creep  strength  or  other  criteria  except  their 
a-c  magnetic  permeability.  The  reduction  of  pole  face  losses 
which  is  offered  by  a  material  which  has  been  selected  for  its 
high  magnetic  permeability  may  eliminate  a  critical  need  for  a 
high  temperature  material  and  the  need  for  special  and  com¬ 
plicated  cooling  passages  in  the  rotor. 

RECOMMENDATIONS 

1.  The  only  ways  by  which  tooth  ripple  losses  can  be 
reduced  with  the  present  design  concepts  are  by  (a) 
increasing  the  radial  gap,  (b)  minimizing  the  amplitude 
of  the  flux  density  variation  on  the  surface  of  the  pole 
face  by  maintaining  the  optimum  ratio  of  slot  width  and 
slot  pitch,  and  (c)  increasing  the  length  or  the  effective 
resistance  of  the  eddy-current  path  by  incorporating 
grooves  in  pole  faces  in  a  direction  parallel  to  the 
motion  of  the  rotor.  The  main  disadvantages  incurred 
by  using  these  methods  for  reducing  the  pole  face  losses 
are  (a)  additional  excitation  is  required,  (b)  leakage 
reactance  in  the  armature  is  increased,  and  (c)  slotting 
rite  pole  faces  is  a  very  expensive  operation  and  only 
increases  the  effective  resistance  of  the  material  by 
a  factor  of  approximately  2.6.  It  is  therefore  recom¬ 
mended  that  other  techniques  for  the  reduction  of  pole 
face  losses  be  investigated.  The  most  promising  solution 
of  the  dilemma  would  be  to  combine  a  material  which  offers 
excellent  magnetic  properties  for  the  reduction  of  eddy- 
current  losses  with  a  material  that  offers  excellent  creep 
strength  and  d—c  magnetic  properties  for  a  solid  rotor. 

In  this  case  the  pole  tips  would  be  made  of  a  laminated 
material  ox  excellent  a-c  magnetic  properties  and  attached 
to  the  pole  bodies.  The  use  of  a  material  that  is  selected 
for  its  high  a-c  magnetic  permeability  about  the  operating 
flux  density  can  reduce  pole  face  eddy-current  losses  by  as 
much  as  7-fold.  This  method  also  permits  a  laminated 
construction  with  interlaminar  insulation  which  provides 
a  large  increase  in  the  effective  electrical  resistance 
to  eddy-currents.  It  is  estimated  that  the  laminated 
construction  can  cause  a  5  to  6 -fold  increase  of  the 
resistance  when  compared  with  the  resistance  of  grooved 
pole  faces.  Thus,  the  eddy-current  losses  can  be  decreased 
by  approximately  2-folds  by  using  a  laminated  construction 
in  place  of  grooved  pole  faces.  The  overall  decrease  of 
eddy-currant  losses  with  a  laminated  material  which  has 
been  chosen  for  its  high  a-c  magnetic  permeability  can 
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approach  16-folds  when  compared  with  losses  associated  with 
material  that  has  been  selected  for  its  creep  strength  or 
other  criteria  except  its  a-c  magnetic  permeability.  The 
greatly  decreased  eddy-current  losses  will  permit  much 
cooler  operation  of  the  rotor  which  will  enhance  the 
strength  of  the  rotor.  The  reduction  in  heating  may  also 
allow  the  rotor  coolant  temperature  to  be  increased  and  even 
eliminate  the  need  for  special  and  complicated  cooling 
passages  in  the  rotor. 

2.  An  estimate  of  the  eddy-current  losses  in  pole  faces 
from  theoretical  considerations  requires  a  knowledge  of 
the  incremental  and  reversible  permeabilities  of  the 
material  in  the  pole.  Therefore,  experiments  should  be 
performed  to  determine  these  properties  for  the  materials 
which  are  candidates  for  a  solid  rotor  and  for  pole  tips. 
Since  H-ll  steel  will  be  used  in  the  experimental,  467  kva 
generator,  it  is  relatively  important  that  the  incremental 
permeabilities  of  this  material  be  evaluated  so  that  a 
comparison  of  theoretical  data  and  experimental  data  can 
be  made.  This  comparison  can  be  used  to  confirm  the 
general  validity  of  the  theoretical  equations. 

The  incremental  permeabilities  of  materials  would  be  evalu¬ 
ated  by  applying  a  steady  induction  on  a  specimen  and  by  apply¬ 
ing  an  alternating  induction  first  in  the  opposite  direction. 

The  incremental  permeability  would  be  determined  from  the  minor 
hysteresis  loops  produced  by  the  alternating  induction.  The 
alternating  induction  would  be  varied  from  zero  to  15%  of  the 
steady  induction.  The  reversible  permeability  would  be  deter¬ 
mined  by  extrapolating  AB  =  0.  The  permeabilities  would  be 
evaluated  along  the  normal  magnetization  curve  and  with  descend¬ 
ing  and  ascending  value  of  the  steady  induction. 

Generally  a  material  which  has  a  high  normal  permeability 
at  a  particular  value  of  induction  also  has  high  incremental 
permeabilities  at  this  level  of  induction.  The  normal  perme¬ 
ability  of  the  cobalt-iron  alloys  at  inductions  up  to  10  kilo- 
gausses  are  relatively  high.  Thus,  these  materials  are  poten¬ 
tially  useful  for  the  reduction  of  pole  face  losses  since  their 
incremental  permeabilities  are  most  likely  high.  It  is  recom¬ 
mended  that  incremental  permeabilities  of  these  materials  be 
evaluated  to  determine  their  capabilities  as  materials  for  the 
reduction  of  pole  lace  losses. 

The  material  which  offers  the  best  combination  of  a-c  mag¬ 
netic  and  strength  properties  for  use  as  pole  tips  should  then 
be  tested  in  the  experimental  generator  for  comparison  with  a 
material  that  has  been  selected  for  the  best  combination  of 
creep  strength  and  d-c  magnetic  properties  for  the  highly  stress¬ 
ed  main  rotor  body. 
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attempt  was  made  to  control  the  temperature  of  the  pole  faces.  There¬ 
fore,  the  test  results  had  to  be  corrected  for  a  change  of  resistivity 
since  the  pole  faces  will  operate  at  a  much  higher  temperature  in  the 
generator . 


A  complete  analysis  of  the  pole  face  loss  problem  based  upon 
theoretical  considerations  is  presented  in  the  report.  The  problem  ofj 
dealing  quantitatively  with  the  magnetic  flux  density  distribution  at  I 
the  surface  of  the  poles  is  treated  analytically.  Harmonic  frequencies 
and  the  losses  due  to  each  harmonic  frequency  were  considered.  The  j 
effect  of  grooves  in  the  pole  faces  were  also  considered  in  the  theore 
tical  analysis.  The  theoretical  equations  show  that  the  eddy-current 
losses  are  a  function  of  permeability  of  the  pole  face  material.  In 
the  actual  problem  this  quantity  is  a  variable  and  a  method  for 
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